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ABSTRACT 

The spectral distribution of field white dwarfs shows a feature called the "non-DA 
gap". As defined by Bergeron et al., this is a temperature range (SlOOiiT-GlOOX) where 
relatively few non-DA stars are found, even though such stars are abundant on either 
side of the gap. It is usually viewed as an indication that a significant fraction of white 
dwarfs switch their atmospheric compositions back and forth between hydrogen-rich and 
helium-rich as they cool. In this paper, we present a Monte Carlo model of the Galactic 
disk white dwarf population, based on the spectral evolution model of Chen and Hansen. 
We find that the non-DA gap emerges naturally, even though our model only allows 
white dwarf atmospheres to evolve monotonically from hydrogen-rich to helium-rich 
through convective mixing. We conclude by discussing the effects of convective mixing 
on the white dwarf luminosity function and the use thereof for Cosmochronology. 

Subject headings: Calaxy: evolution — white dwarfs — methods: numerical 



1. The non-DA gap as a complication of WD Cosmochronology 



White dwarfs (WDs) are the endpoint of stellar evolution for low mass stars. Their spectra can 
be divided into two broad categories: the DA stars (~ 75%), whose atmospheric composition is ob- 
served to be hydrogen- dominated, and the non-DA stars (~ 25%), whose atmospheric composition 
is helium-dominated. The evolution of WDs is mainly driven by the cooling of the degenerate core. 
As a result, the effective temperature (Tg//), or equivalently, the luminosity of a W D can serve as a 
cosmic clock, indicating the time passed since its format ion (IFontaine et al.ll2001h. Analytical an d 



numerical models of WD cooling have been made (e.g., 


Mestel ( 


1952 




Lamb & van Horn ( 


1975): 


Wood ( 


1992 


): 


Althaus & Benvenuto ( 


1997): 


Hansen (19991): 


Chabrier et al. ( 


2nnnl): Salarls et al. 


(2nnd): 


Garcfa-Berro et al.l ( 


20081^: 


Renedo et al. 


( 


20ld)), allowing the use of white dwarf luminos- 



ity functiori (WDLF) to derive the age of Galactic disk (Liebert et al. 1988 : [winget et al 



2007 



Wood I1992I) and other stellar populations ([Hansen et al. 



Bedin et al.ll2009l ). The practice is known as WD Cosmochronology (jWinget et al.lll987l ). 
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However, there are still many coni p lications th at impede WD Cosmochronology. Among them 



are the observations of Bergeron et alj (|l997l . l200ll ). In their work, photometric and spectroscopic 



observations of 152 cool (4000i^ < Teff < 12000i^) WDs that belong to the local Galactic disk 
are analyzed with state-of-the-art model atmospheres of pure hydrogen, pure helium and H/He 
mixtures. It was shown that the ratio of observed helium-rich stars to hydrogen-rich stars is 
greatly reduced in the temperature range 5100if-6100iC and enhanced both above and below this 
temperature range (see Figure [T]) . This temperature range was consequently termed the "non-DA 
gap" to indicate that few non-DA stars are found within. Furthermore, some so-called "peculiar 
stars", whose spectral distributions are better reproduced with hydrogen model atmospheres but 
which show no Ha lines, are found in the temperature range ~ 6000A' to ~ 8500X. The results 
seemed to portray the following picture of WD spectral evolution as it cools from Tgjj = 12000-fC 
to Tgff = 4000i^: First, a significant fraction of the DA stars (which are the majority of all WDs) 
switch their atmospheric composition from hydrogen-rich to helium-rich in the temperature range 
12000K ~ SSOOK, causing a steady increase of non-DA ratio (with respect to the decrease of T^ff). 
Later, some helium stars start to evolve into peculiar stars (at ~ SSOOiT), and then into hydrogen- 
rich stars (at ~ GlOOiiT). The stars remain hydrogen-rich in the "non-DA gap" (GIOOX ~ 5100-R') 
and transform into helium-rich WDs again at ~ SlOO-fC. The observations introduced an uncertainty 
in WD Cosmochro nology because atmospheric composition significantly affects cooling rate at lower 
temperatures (e.g., iHansenl (j 19991 )). yet, all of the previous works on theoretical WDLF were done 
with only models of fixed atmospheric composition. 

Bergeron et al. ( 1997 ) proposed an accretion-convective mixing model to account for the 



observations. In their model, hydrogen-rich inter-stellar medium is first accreted onto the surface 
of a helium WD at temperatures higher than 6100-fC without being spectroscopically visible, due to 
the quenching effects of the high photospheric pressure that characterizes these atmospheres. Such 
a WD will appear as a peculiar star, until the hydrogen concentration reaches the level at which 
starts to recombine. At this point, a runaway process will occur and turn the star into a DA star 
at ~ GlOOiT. The spectral type of the star remains DA as it cools from GlOOiT to SlOOiT. Finally, 
at SlOOiT, it is transformed into a helium-rich star again through convective mixing. However, 
this proposal was later rejected bv iMalo et al.l (|l999l ) as they claim no models in their grids were 
found to have Ha lines quenched while retaining the spectral energy distribution of a pure hydrogen 
atmosphere. 



An alternative explanation was given in lHansenI (| 19991 ). where the non-DA gap was suggested to 
be the consequence of t he different cooling rat es of hydrogen-rich and helium-rich WDs. However, it 
was later pointed out in [Bergeron et al.l (|200ll ) that such a proposal cannot explain the enhancement 
of the non-DA ratio at temperatures above the blue edge of the gap. As a consequence, there is 
still no convincing explanation for the non-DA gap and its properties. 



^We refer the reader to section [2] for a review of convective mixing. 
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2. Convective Mixing as a Mechanism of Spectral Evolution 



It has been estabhshed that elemental diffusion is efficient in the WD atmosphere (jMuchmore 



19841 ) ■ In the absence of other mechanisms, a WD with only small amounts of hydrogen therefore 
possesses a photosphere close to pure hydrogen. The leading mechanism to transform a hydrogen- 
rich (photos phere) WD to a helium-rich WD is through a process called convective mixing, first 
proposed by IShipmanI (| 19721 ) . The idea is as follows: The surface convection zone of a hydrogen 
WD grows as T^^jf decreases. For WDs with intermediate amount of hydrogen (i.e., the hydrogen 
layer is optically thick but less massive than W~^Mq), the convection zone will eventually reach 
the underlying helium layer. The helium will thus be brought to the surface through convective 
motion, changing the photospheric composition. 

Th is mechan i sm ha s bee n investigated extensively , and a number of outcomes have been pro- 
posed. IShipmanI (jl972l ) and iBaglin fc Vauclaiii (jl973l ) suggested that a temper ature increase ac- 
companies convective mixing. However, the evolution model of iKoesteii (|l976l ) appears to show 
no such behavior. Besides, while most groups conjectured (based on the size of convection zones 
in pure hydrogen and helium WDs) that th e composition of the resultant photosphere would be 
close to pure helium, [Bergeron et al.l (j 19901 ). based on observational results, suggested that the 
convective m ixing might resul t in an H/He mixture whose opacity is still dominated by hydrogen. 
The work of I Chen &: HansenI (j201ll ) succeeded in unifying all of the suggested outcome into one 
framework. Using WD mass (Mwd) and total hydrogen mass {mn) as the input parameters, a set 
of self-consistent, quantitative cooling curves and chemical evolution curves were derived (see, e.g.. 
Figure [2]) . It was concluded that the temperature increase which accompanies convective mixing 
only occurs when the post-mixing envelope is convectively coupled to the degenerate core, and 
that the post-mixing photosphere hydrogen fraction {Xgurf) is a function of both M\yd and ruH- 
It has been shown that massive WDs with a higher amount of hydrogen are more likely to stay 
hydrogen-rich after convective mixing. 

The conve c tive m ixing mechanism has also been deployed in the aforementioned model of 
Bergeron et al.l (| 19971 ). in which the authors seemed to implicitly assume no temperature increase 
and a post-mixing ph otosphere which can be treated as pure helium. Based on the results of 
Chen fc HansenI (120111 ). these assumptions require some modifications. The T^jf should increase 
upon mixing because the WD envelopes are already convectively coupled to the degenerate core. 
Another feature of the convective mixing WDs that Bergeron et al. did not fully exploit is their 
faster, post-mixing cooling rate (cf. Figure [2]). 

In the rest of this paper, we will demonstrate that the non-DA gap naturally emerges from 
the improved treatment of convective mixing, without additional mechanisms. In order to do this, 
we will construct a model realization of the Galactic disk WD population based on the recent 
determinations of Galactic parameters. 
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3. A Monte Carlo realization of the Galactic white dwarf population 



In this section, we will describe the construction of a Monte Car l o sani ple of Galactic WDs 
which includes the treatment of convective mixing of IChen &: HansenI (1201 ih. The purpose is to 
see if such inclusion could account for the non-DA gap of Bergeron et alj (jl997l . l200lh . Specifically, 
the non-DA to total WD ratio as a function of T^ff will be calculated at a specific Galactic disk 
age td- We will also compare the WDLF synthesized from this sample to some other theoretical 
and observational WDLF and discuss the cosmological implications. 

The Monte Carlo sample is defined by an array (~ 10® entries, one entry per star) that specifies 
the physical properties of each WD (or its progenitor, if the star is too young to have evolved). Two 
random numbers, stellar formation time (tb) and main sequence stellar mass (M^s), are drawn for 
each star. As an approximation, we take the star formation rate of our Galaxy to be constant and 
generate tf, to be unifo rmly distribute d between and td- Mms is generated to conform with the 
initial mass function of iKroupal (120011 ) with an upper limit of 8Mq (stars with masses > 8Mq are 
not statistically important). A pr e- WD lifespan t Mfi is assigned to each star based on its Mms-, 



following the formula described in iHurley et al.l (2000). A comparison between td and {tb + tMs) 
is then made: If td < (tf, + tjus)) the star is not yet a WD and is discard ed; If td > {tb + tMs)i '^^ 
further calculate M\yd using the relation established bv lZhao et al.l (j2012l ). We define the difference 
between the two quantities as the cooling age tc, i.e., tc = td — tb — tMS- 

All physical quantities of a WD can be derived once Mw^, tc, and a WD evolution model 
grid are specified. Three WD evolution grids are available in this work: the pure helium WD 
models with mn being equal to zero, the thin hydrogen layer WD models with tuh = 1O~^M0 
and the thick hydrogen layer WD models with a hydrogen fractional mass q{H) = mn /Mw£) = 
10~^. All WDs in our sample have helium fractional mass q{He) = 10~^. The pure helium WD 



models and thick hydrogen layer W D models are des c ribed in lHansenl (|l999l ) and the thin hydrogen 
layer WD models are described in IChen &: HansenI (j201ll ). The mn value we chose for the thin 
hydrogen layer stars, IO^^Mq, is a value consistent with asteroseismology studies. The hydrogen 
fractional mass we chose for the thick hydrogen layer stars and the helium fractional mass we chose 
for all WDs are the maximum amount that ca n survive the pre- WD evo l utionary phases and is 



Fontaine et al 



2001 



considered as the "standar d thick layer model" (iMichaud &: Fontaine 
1990 



1984 



DAntona fc Mazzitelli 



). It is a simplification to account for all WDs with only the mentioned 
three model grids. In any real WD samples, either mn or q{H) should more or less follow a 
continuous distribution. The mn or qh values we adopt in our work can be viewed as a set of 
characteristic values, i.e., where the distribution is peaked. 

In the following paragraphs, we show the results of one of our simulations with td = 9.5 Gyr. 
We assume 15% of the WDs are pure helium WDs, 32% of the WDs are thin hydrogen layer 
WDs and 53% of the WDs are thic k hydrogen layer WDs. The choice of disk age, td = 9.5 Gyr, is 
consistent with lLeggett et al.l (jl998l ). The WD type fractions are chosen to match the spectral type 
distribution from observations (cf. Figure H]). However, it should be kept in mind that our main 
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purpose is to show the quahtative consistency between our model and observations. No attempts 
were made to accurately fit for either the Galactic age or the percentage of convective mixing WDs. 

Figure [3] shows the distribution of a representative Monte Carlo WD population in the T^ff- 
MwD space and Tefj-Xsurf space. We have reduced the number of stars in the sample to avoid 
the plot being over-congested (we used the full amount of numerical stars in making Figure H] 
and Figur e [3 and con f irmed that the results have already converged with respect to sample size). 
Following iMalo et al.l (jl999l ). WDs with Xgurf < 10~^ are classified as helium-rich stars and those 
with Xsurf > IQ-^M© are classified as hydrogen-rich stars. Interestingly, a small population of 
intermediate Xg^rf stars (10~^ < X^urf < 10 ~^) are found to cluster arou nd 8500i^ to 6000K, 
which is the location of the peculiar stars of lBergeron et al.l (jl997l . l200ll ). Therefore, they seem to 
be natural candidates for the observed population of peculiar stars. These stars are evolved from 
intermediate massive (O.69M0 < MwD ^ 0.87 Mq) thin hydrogen layer WDs. The higher stellar 
masses of these WDs make their convection zones shallower (compare to the majority of WDs which 
are less massive) and consequently leads to the intermediate hydrogen fraction that results upon 
convective mixing (those with even higher masses would remain hydrogen-rich stars). At later 
times, most of them evolve into true helium-rich stars due to the deepening of their convection 
zones, hence they are less likely to be found below GOOO-RT. We would like to emphasize that ou r 
plot (upper panel, the T^jf-My/D plot) is not a simulation of Figure 21 of lBergeron et al.l (j200ll ) 
(i.e., our Figured]) and is intrinsically different from it in the following two respects. First, the 
"stellar density distribution" is different (the stars in Figure [3] are aggregated in the lower-right 
region whereas the stars in Figured] are more uniformly distributed), because our plot is made with 
the whole sample instead of the detected sample. In order to properly compare observational results 
with our numerical work, information on detectability needs to be incorporated (cf. the next two 
paragraphs). Secondly, our numerical sample does not have any WDs less massive than ~ O.SSMq, 
because the progenitor of these low mass WDs would still be on the main sequence track, according 
to our Mms~Mwd relation. The low mass WDs in Figure [1] could be evolved from, e.g., binary 
systems, and are not accounted for in our model. 



While the non-DA gap is qualitatively apparent in iBergeron et al.l (120011 ) . the most up-to- 
date and quantitative obser vation al data on th e non- DA ratio of cold WDs is in the works of 
Tremblay k, Bergeron! (j2008l ) and iLeggett et al.l (jl998l ) . In Figure IH we will dir ectly compare 
our numerical results with these observations. The data in IBergeron et al.l (120011) ha ve a larger 
sample size, however, detectability information, e.g., the Vmax correction ( Schmidtl 1975 ). is absent. 
Therefore, a quantitative comparison with our work cannot be made at this point. 

Two sets of observational data are plotted in Figure |4] in comparison with our numerical 
results (solid line). The two dotte d lines are the upper and low er limits of the non-DA to total 



WD ratio derived from Figure 11 of lTremblav &: BergeronI (20081 ). The filled circles with error bars 



are compiled by us from the results of lLeggett et al.l (| 19981 ): instead of b inning their WD sam ple 
(which is an updated determination of the 43 faintest WDs presented in iLiebert et al.l ( 1988 )) in 
luminosity, we binned the sample (with Vmax correction) in Tg j f . The binning was done for the DA 
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(hydrogen-rich) stars and the non-DA (helium-rich) stars respectively. In each bin, we calculate 
the upper limit of the non-DA to total WD ratio by dividing the upper limit of the non-DA density 
to the lower limit of total WD density (while fixing the non-DA density at its upper limit), and the 
lower limit of the ratio by dividing the lower limit of the non-DA density to the upper limit of the 
total WD density (while fixing the non-DA density at its lower limit). All of the results are plotted 
in Figure m except that corresponds to the bin SOOOX-QOOO/C, where the result and statistical error 
are of the same size (because there is only one star in the sample). Our numerical result (solid line) 
is generated by directly dividing the number of non-DAs to the number of all WDs in each T^ff 
bin (lOOOK in size) of our Monte Carlo sample. We can see that a non-DA gap naturally emerges 
in the range SOOO-ftT-GOOOiT as a consequence of convective mixing and the resulting faster cooling 
rates of post-mixing WDs. 

In Figure [U we examine how this qualitative change in the evolution of individual stars is 
reflected in the luminosity function of the population as a whole. We compile the WDLF from 
our Monte Carlo sample which includes convective mixing WDs (CM- WDLF hereafter, solid line), 
together with that from a sample composed solely of thick hydrogen layer WDs (TH-WDLF here- 
after, dash-dotted li ne) and comp a re to three different sources of observation al WDLF. Th e three 



different sources are iLiebert et al.l (|l988l ) (LDM hereafter, for L > IQ-'^-^L.t.). iLeggett et all (119981) 



(LRB hereafter, for L < 10 Lq), and the results from Sloan Digital Sky Survey ([Harris et al 



20061 ) (SDSS hereafter, for the entire range of luminosity) . 



The CM- WDLF starts to deviate from the TH-WDLF below L = W^^Lq. The value of the 
CM- WDLF is larger than that of the TH-WDLF between L = lO'^L© and L = W'^-^Lq, due 
to the temperature increase that accompanies convective mixing. Below L = lO^'^'^L©, however, 
the situation inverts, because the convective mixing WDs (and the pure helium WDs) have faster 
cooling rates, resulting in a flatter, less peaked WDLF. The faster cooling rates also enable the 
WDLF to extend to a fainter luminosity. The observations from LDM and LRB are not sufficiently 
densely sampled to distinguish between the two different theoretical shapes, but the SDSS sample 
may have that capability once its properties are more accurately determined (e.g., when the parallax 
of each star are obtained to determine the distances; currently, the SDSS WDLF are derived with 
photometric distance and an idealized assumption of logg = 8.0 for all WDs is made). 



4. Conclusions and discussions 



We have demonstrated that the non-DA gap can be reproduced i i i a p opulation model by 
exploiting the spectral evolution model investigated in I Chen &: Hansen! (j201ll ) . The non-DA gap 
can be formed by the combined effects of convective mixing and the higher cooling rates of the post- 
m ixing WDs. An evolut ion model that transforms DC stars into DA stars, e.g., the one outlined 



m 



Bergeron et al.l (jl997l ). is not needed to produce the non-DA gap. 



Spectral evolution is important to WD Cosmochronology because the cooling rate of a WD is 
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significantly affected by tlie atmosplieric conditions (cf . Figure [2]) . Tlie nature of the faint tail of 
the observational WDLF based on the SDSS sample could be indicative of the problems faced by 
theoretical WDLF made with only thick hydrogen layer WDs. 

Convective mixing, the non-DA gap, and WD Cosmochronology are closely related, in that 
advances in each topic will bring new insights to the others. While this study has demonstrated the 
relation between them and established the consistency within the available observational results, 
works on each front, both theoretically and observationally, are not complete. To that end, a 
spectroscopically confirmed WD sample larger than that of LDM/LRB is in need. Although the 
sample from SDSS represents an important advancement, it suffers from the selection effects and 
halo WD contamination. The results from future WD surveys will tell us more about the history 
of our Galaxy. 

The authors acknowledge D. E. Winget, M. H. Montgomery, S. O. Kepler, and J. Adamczak 
for discussion and helpful comments. 
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Fig. 1. — A simplified version of Figure 21 of iBergeron et alJ ()200ll ). Filled Circle: H-ricli star, 
Hollow Circle: He-rich star, Triangle: peculiar star. The main atmospheric composition is taken 
from their Table 2. Complications such as DQ stars, C2H stars. . . , etc. are not presented. 
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Fig. 2. — An e x ampl e of the cooling curves and chemical evolution curves of 0.6Mq WDs in 



Chen fc HansenI (|201lh (Figures 11 and 15, © 2011 RAS). UPPER: The cooling curve of the 
convective mixing WD with mn = 1O~^M0 (solid line), that of the thick hydrogen layer WD 
(dashed line), and that of the pure helium WD (dash-dotted line). LOWER: Chemical evolution 
curve of the convective mixing WD. In the other two WDs, no chemical evolution takes place. 
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Fig. 3. — Distribution of our representative Monte Carlo WD sample in two different parameter 
spaces. The use of symbols is identical to that in Figure [TJ UPPER: Distribution in the T^fj-MwD 
space. The dash-dotted line marks the T^ff where convective mixing first occurs. The dashed line 
marks the locus where the newly-mixed WDs converge. The two dotted lines bracket a region of 
low non-DA ratio, corresponding to the "non-DA gap". We note the T^fj evolution of WD is non- 
monotonic when convective mixing occurs. LOWER: Distribution in T(,f f-Xg^rf space. The dashed 
line marks the locus where the WDs converge after convective mixing. The locus where convective 
mixing first occurs coincides with the upper edge of the plot and is therefore not explicitly present. 
The thick hydrogen layer DA stars remain on the upper edge throughout their lifetime. The pure 
helium WDs do not appear as they are located at infinity on this log scale plot. The WDs cool 
faster after convective mixing and are able to reach the low Tg j j end of the plot in large quantities 
within the given Galactic age. 
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Fig. 4. — Comparison of the non -DA to total WD ratio from our numerical work (solid lin e ) to t he 



observational data (dotted lines: iTremblav Sz Bergeron 
The binning intervals we used to compile the data in 

[5000/C,6000i^]. . . , etc., to conform with the data of Eremblav BergeronI (|2008l ) 



20081^. filled ci rcles: 



Leggett et al.1 (jl998|)) 



Leggett et al.l (119981 ) are [AOOOK, 5000/C] 
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Fig. 5.— UPPER: Theoretical CM-WDLF (solid lin e) and TH-WDL F (dash-dotted line), together 
with observational WDLF from lLiebert et aL ( 1988 ) (filled circle) and Leggett et al. ( 1998 ) (hollow 
circle). LOWER: Same theoretical construction compared with observational WDLF from SDSS 
([Harris et ahlbood ). 



